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ABSTRACT 

The Australia Telescope Compact Array (ATCA) has been used to image class I methanol 
masers at 9.9, 25 (a series from J=2 to J=9), 84, 95 and 104 GHz located in the vicinity of 
IRAS 16547-4247 (G343. 12-0.06), a luminous young stellar object known to harbour a radio 
jet. The detected maser emission consists of a cluster of 6 spots spread over an area of 30". 
Five spots were detected in only the 84- and 95-GHz transitions (for two spots the 84-GHz 
detection is marginal), while the sixth spot shows activity in all 12 observed transitions. We 
report the first interferometric observations of the rare 9.9- and 104-GHz masers. It is shown 
that the spectra contain a very narrow spike (< 0.03 km s~^) and the brightness temperature 
in these two transitions exceeds 5.3 x 10^ and 2.0 x lO"* K, respectively. The three most 
southern maser spots show a clear association with the shocked gas traced by the H2 2.12 jim 
emission associated with the radio jet and their velocities are close to that of the molecular 
core within which the jet is embedded. This fact supports the idea that the class I masers reside 
in the interface regions of outflows. Comparison with OH masers and infrared data reveals 
a potential discrepancy in the expected evolutionary state. The presence of the OH masers 
usually means that the source is evolved, but the infrared data suggest otherwise. The lack of 
any class II methanol maser emission at 6.7 GHz in the source raises an additional question, 
is this source too young or too old to have a 6.7 GHz maser? We argue that both cases are 
possible and suggest that the evolutionary stage where the class I masers are active, may last 
longer and start earlier than when the class II masers are active. However, it is currently not 
possible to reveal the exact evolutionary status of IRAS 16547-4247. 
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1 INTRODUCTION 

Methanol masers are commonly found in massiv e star-forming re 



gions. They fall into 2 categories first de fined by 
il988l) and further considered by Voron kov et al. 



Batrla & MentenI 
j2005bh : clasTn 



masers (of which the 6.7 GHz is the best known and usually 
strongest) are closely associated with infrared sources and reside in 
the close environment of exciting (proto-)stars. In contrast, class I 
masers (e.g. 95 GHz) are usually found offset (by up to a par- 
sec) from continuum sources (e.g., Kurtz, Hofner & Alvarez 2004). 
Theoretical calculations strongly suggest that class I masers are 
pumped via collisions with the molecular hydrogen i n contrast to 
class II masers, which have a radiative pum p (e.g., ICragg et'all 
1 1992tlVoronkovll 99911 Voronkov et alj2005d) . 

Class I masers are relatively poorly studied. There is direct 
observational evidence to suggest a relationship between these 
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masers and outflows in a number of sources jPlambeck & MentenI 
ll990';'Menten"l99^, although another scenario, which in volves 
cloud-cloud collisions , has also been proposed lSobole\^ ll992t 
iMehringer & Me nten 1995 ISaUi. Sobolev & Kalinina 2002). The 
common point of these two scenarios is the presence of shocks. 
There is observational evidence that methanol abund ance is signifi- 
cantly increased in the shock processed regions (e.g. jGibb & Davis! 
1998) and the gas is expected to be heated and compressed in such 
regions providing more frequent collisions and, therefore, more ef- 
ficient pumping. 

To further investigate the relationship of different class I maser 
transitions with outflows and to provide data for theoretical stud- 
ies, we have made interferometric observations of 12 class I maser 
transitions towards IRAS 16547-4247 (G343. 12-0.06) using the 
Australia Telescope Compa ct Array (ATCA). This is a sign ificant 
fraction (50% according to lMuller. Menten & MadeJl2004ll) of all 
known class I maser transitions, the majority of which belong to 
the J2— Ii E series. We observed at least one transition from each 
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Table 1. Dates of observations, array configurations and observed transitions 
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transi tion series known to exhibit class I maser activity ^Voronko^^ 
Il999h and covered all known class I maser transitions within the 
frequency range of existing ATCA receivers. Only the J=5 tran- 
sition from the J2— Ji E series at 25 GHz has been pr eviously 
observed towards IRAS 16547-4247 iVoronkov et aljEoO Sa). al- 
though with no position measurement. In this paper we present the 
first detection of 7 other transitions from this series as well as the 
first detection of the 84-GHz maser towards this source. The other 
3 transitions were detected in a number of single dish and ATCA 
maser line surveys (Slvsh et al .il994lw^tts et all200d Voronkov 
et al. unpublished observations). It is worth mentioning that this 
source is the only example where all these class I maser transitions 
have been detected together, including the relatively rare 9.9 and 
104 GHz masers. 

IRAS 16547-4257 has recently been the focus of attention as 
the most luminous young stellar object (YSO) known to harbour a 
radio jet. At a distance of 2.9 kpc (Bronfman, private communica- 
tion), the YSO has a bolometric luminosity of 6.2 x 10* L©, equiv- 
alent to that of a single 08 Z ero Age Main Seque nce (ZAMS) star. 
The jet was first detected bv lGarav et alJ j2003|) using the ATCA 
at 1.4, 2.5, 4.8 and 8.6 GHz and subsequently observed using the 
VLA at 8.5 GHz and 14.9 GHz ( Rodn'auez et al. 2005). This is 
the first reported case of a radio jet associated with a young O- 
type star. Located at the centre of a massive dense molecular core 
(approximately 10^ M0), the jet is driving a hi ghly energetic colli- 
mated bipolar outflow that extends over 1.5 pc iBrooks et alj2003l : 
barav et al,2006.) . 



2 OBSERVATIONS 

Observations were made using the ATCA in 2005 May-August 
with assorted array configurations and various correlator modes. 
The details, including the spectral resolutions, velocity coverages, 
full-width at half maximum (FWHM) of the primary beam, and 
synthesized beam parameters, are summarized in Table Q We 
adopted the rest frequencies determined by Miiller et al. ( 2004). 
The measurement uncertainties are given in Table Q in parenthe- 
ses after each corresponding frequency. They are expressed in units 
of the least significant figure. Corresponding uncertainties in the 
radial velocity are given in the next column. 



The J=2,3 and 4 transitions of the J2— Ji E series near 25 GHz 
were observed in a single 8 MHz band centred on the average fre- 
quency of these transitions. The remaining transitions were ob- 
served using a separate frequency setup for each individual tran- 
sition. The correlator setup used for observations of the J=5 and 6 
transitions on May 2 allowed us to observe the 24896 MHz con- 
tinuum emission with a 128 MHz bandwidth simultaneously with 
the spectral line observations. However, for the subsequent obser- 
vations we chose a higher spectral resolution rather than a con- 
tinuum channel. With the exception of the 9.9 GHz observations 
on June 16 two orthogonal linear polarizations were recorded. For 
June 16 only one linear polarization was recorded in order to obtain 
a larger number of spectral channels and achieve a spectral resolu- 
tion comparable with that of the higher frequency t ransiti ons. It is 
worth noting that lWiesemever. Thum & Walmsle^l i2004) found a 
number of mm-wavelength class I methanol maser transitions to be 
polarized. The correlator setup used in our study did not allow us to 
determine polarization properties. However, during the data reduc- 
tion both polarizations, if present, were averaged together. Hence, 
all the flux densities with the exception of that at 9.9 GHz measured 
on June 16 are not affected by polarization. For the latter we exam- 
ined the flux density measured from 10 individual scans taken at 
various hour angles assuming that the source was unresolved and 
located at the position found by the imaging. No parallactic an- 
gle dependence was found beyond the noise level of 0.3 Jy (Icr). 
This fact as well as a comparison with the May observations (see 
section 1331 performed in the dual polarization mode leads us to 
conclude that polarization effects, if present, are less than the cali- 
bration uncertainty. 

The position of the phase and pointing centre was 02000 = 
16''58™16';57, (52000 = -42°52'24'.'0 in all observing sessions. 
Reference pointing using the source 1646-50 (which served also as 
a phase calibrator), was used for all observations except those at 
9.9 GHz, where a global pointing model was used. From the statis- 
tics of the pointing solution the reference pointing accuracy was 
estimated to be 2'.'3±(/.'8 for all observing runs. The accuracy of 
the global model (which affects the 9.9-GHz data only), was esti- 
mated to be 6"±3" and 10" ±4" for the 2 May and 16 June runs 
respectively. The accuracy of the reference pointing affects the ac- 
curacy of flux density measurements, particularly for sources offset 
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Figure 1. The general morphology of the region: the 95 GHz maser spots and the 25 GHz continuum sources are overlayed on the H2 2.12/im image 
tBrooks et al. 2003) . The continuum sources are also shown on the large scale image and believed to delineate the direction of the jet. An artifact in the 
northern part of the large scale H2 image is caused by a star. Contours in the small scale image are (1, 2, 3, 4, 5, 7, 9) x3.6 Jy beam^^ and (0.3, 0.5, 1, 3, 5, 7, 
9) X 1.5 mjy beam"! for the 95-GHz maser and 25-GHz continuum sources respectively. The lowest contour of the 25-GHz continuum emission corresponds 
to la. The noise level of the 95-GHz image varies significantly across the image (see Tablelsl. For spot F the lowest contour corresponds to 9a. The numbers 
in parentheses are approximate peak radial velocities in km s~^ for each maser spot. The synthesized beam is shown in the bottom-left comer of the inset 
(small and large ellipses for 25 GHz and 95 GHz, respectively). The open box and two crosses near the central continuum source mark the positions of the 
mainline OH masers and two well separated clusters of the H2O masers, respectively. 



from the pointing centre. The positional accuracy of the maser loca- 
tions is primarily influenced by the quality of the phase calibration 
and is believed to be better than 0.5 arcsec. All 3mm observations 
were undertaken as two-point mosaics. An additional pointing was 
offset half of the primary beam width (see Table0 to the East from 
the position given above. 

Data reduction was performed using the MIRIAD package (28 
August 2005 release) following standard procedures with the ex- 
ception of bandpass calibration. Using the UVLIN task of MIRIAD a 
low-order polynomial fit to the spectrum of the bandpass calibrator 
(1921-293) has been performed before attempting a bandpass so- 
lution with the MFCAL task. This approach enables us to achieve a 
higher spectral line dynamic range (i.e. detect a weaker spectral line 
in front of the continuum of a given flux density), although none of 
the datasets described in the paper turned out to be dynamic range 
limited. 

The flux density scale at 9.9 and 25 GHz was established us- 
ing the standard ATCA primary calibrator, 1934-638. For higher 
frequencies planetary calibration was done using Uranus. In order 
to assess the variability of the source some transitions were ob- 
served more than once. Two such repeated observations (see Ta- 
ble0 were quite short and no imaging was possible. However, an 
accurate check of the flux density scale is still possible provided 
the source position is the same for all epochs. The accuracy of the 
absolute scale of flux density is estimated to be 3% and 10% at 



9.9 and 25 GHz, respectively^. The uncertainty of high frequency 
planetary calibration of the ATCA has yet to be fully investigated 
and may be as large as 50%. However, based on our experience 
we believe it to be better than 30% for the calibration and data re- 
duction procedure we followed in our observations. This figure is 
in agreement with the measured relative flux densities of the sec- 
ondary calibrator at different frequencies. An additional 8% flux 
density uncertainty related to the pointing accuracy given above ex- 
ists for the high frequency (84 and 95 GHz) features situated near 
the edge of the imaged field of view. 



3 RESULTS 

3.1 General morphology 

Figure [T] shows an H2 2.12 fim image of IRAS 16547-4247 
jBrooks et all2003h . The image contains a complex chain of emis- 
sion with three major concentrations consistent with the morpho- 
logical characteristics of Herbig-Haro (HH) objects arising from 
the interaction of a co llimated flow with the ambient medium 
jReipurth & BdivlbOOll) . There are a number of mechanisms re- 
sponsible for excitation of 2.12 fim emission, the most common 

^ For details on the calibration using 1934-638 query this calibrator at the 
ATCA calibrators webpage \http://www.narrabri.atnf.csiro.au/caUbrators'^ 
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of which are fluore scence and coflisional excitation in shocks (e.g. 
iHabart et alJl2005l) . Several mechanisms can simultaneously con- 
tribute to the observed emission. In general, a detailed spectro- 
scopic analysis is required to determ ine the exa ct excitation condi- 
tion (e.g., Eisloffel et al. 1996; Wals h. Lee & B urton 2002). How- 
ever, the morphology of the region, which suggests an interaction 
between the outflow and the ambient medium, together with the 
lack of Br7 emission which indicates that the number of u ltra-violet 
photo ns is insufficient to pump the fluorescence (Brooks et alj 
l2003h . supports the idea that most of the H2 emission in this source 
traces shocked gas. 

The outer concentrations of H2 are located approximately 
symmetrically offset from the brightest source of the cm- 
wavelength radio continuum emission (see Fig.U . This continuum 
source was first detected bv lGarav et alj i2003l> and interpreted as 
thermal (free-free) emission from a radio iet. iGarav et alj i2003h 
also detected two satellite continuum sources and interpreted them 
as the internal working surfaces of th e jet. Sensitive h igh spatial 
resolution observations carried out bv lRodriguez et aP ^2QQ5^ re- 
solved the northern lobe into a chain of spots aligned along the jet 
direction and showed that the peak spot most likely corresponds 
to optically thin thermal emission. The southern lobe was also re- 
solved by Rodriguez et al. (2005), but in contrast to the northern 
lobe a non-thermal origin of the brightest component has been con- 
firmed. Both thermal and non-thermal emission mechanisms are 
expected according to theoretical studies of electron ac celeration 
in shock waves (e.g. lHenriksen. Ptuskin & Mirabel 199 ih . 

In addition to these sources. [Rodrfguez et alj J200^ detected 
three continuum sources displaced by a few arcseconds from the jet 
axis, which is thought to be a line joining the outer continuum lobes 
(Fig. 0. The brightest two of these sources are located between 
the southern non-thermal lobe and the central source and form a 
wiggly structure, convex in the opposite direction from the bow- 
shaped middle H2 concentration. iRodrfguez et alJi200^ proposed 
that one of these continuum sources (the most distant from the jet 
axis) could be associated with a young, low-mass star with a gy- 
rosynchrotron emission. If this interpretation is correct, the curved 
shape of the middle H2 concentration could be the result of an ad- 
ditional outflow from this less massive star. However, it follows 
from the large scale distribution of the H2 emission shown in Fig.Q 
that the molecular outflow covers a rather wide area along the jet. 
Therefore, this morphology could also be caused by a number of 
factors su ch as irregularities in the ambient medium or precession 
of the jet jBrooks et all2003h . 

The expanded scale inset of Fig.Qshows the location of all de- 
tected maser spots at 95 GHz (labelled A to F). The 25 GHz contin- 
uum emission also shown in the figure closely resembl es the triple- 
source morphology found at lower frequencies by iGarav etTH 
i2003h . This continuum measurement and the spectral energy dis- 
tribution (SED) of the source are discussed in a separate paper 
(Brooks et al. 2006). The outermost lobes of the H2 emission are 
outside the ATCA field of view at 3 mm. The southern maser spots 
(labelled A to C) show a remarkable correlation with the promi- 
nent H2 feature (Fig. [T}. Contrary to our e xpectations based on 
the single-dish spectra of lVal'tts et alj i200Cl) obtained at five off- 
set pointings, no maser spots offset in the East-West direction were 
detected. Instead, a group of spots (E, D and F) was found near 
the edge of the field of view in the northern direction. Most likely 
this discrepancy is related to the pointing model used a t the Mopra 
radio telescope during the observations of lVal'tts et alj 



3.2 Methanol masers 

For each spot and observed frequency we produced spectra by sum- 
ming the flux density in each plane of the deconvolved cube after 
the primary beam correction or linear mosaicing. The summation 
was over a box centred on the peak pixel in the image. The size of 
this box was chosen individually for each spot to suppress all unre- 
lated emission (which is difficult at the edge of the field of view), 
and was typically about 4". The spectra for all spots other than 
spot B, are shown in Fig.|2| The horizontal error bars on these plots 
show the 3(T velocity uncertainty related to the uncertainty of the 
rest frequency. With the exception of spot B, all were detected only 
at 84 and 95 GHz, although the 84-GHz detections for spots D and 
E are marginal (both at the level of 4ct). 

The profiles look similar at both frequencies for spot A. How- 
ever, the two spectral features of spot C show significantly differ- 
ent 95-GHz peak flux densities, while at 84 GHz they are similar. 
The 84-GHz profile of spot D looks somewhat similar to its 95- 
GHz counterpart, giving us confidence that this marginal detection 
is real. The case for spot E is more doubtful, a narrow spike near 
—30.33 km s~^ at 84 GHz may be real given the presence of a fea- 
ture with a similar appearence at 95 GHz. The 84-GHz spectrum 
also shows a hint of broad line emission centred at —30.2 km s~^ , 
with a narrow line at —32.18 km s^^ which does not have a 95-GHz 
counterpart. We have formed similar spectra for adjacent areas in 
the image to check that these 84-GHz features are present towards 
the spots D and E only. However, we can not completely exclude 
the possibility that they are spurious because spots D and E are 
located near the edge of the mosaic, and the source contains emis- 
sion at different locations with overlapping velocity ranges. In the 
image cube these features can not be distinguished from the edge 
noise. Therefore, we have not drawn any conclusions from the pos- 
sible presence of 84-GHz emission towards spots D and E. The 95- 
GHz spectrum of spot E contains several components well spread 
in velocity, although the image of the spot E is simple (see Fig.0. 
This is not observed for any of the other spots in this source. The 
presence of multiple components for spot F is more pronounced at 
84 GHz. The 95-GHz profile of this spot looks like a single asym- 
metric broad feature. This may indicate fine scale structure, which 
differs at these two frequencies (i.e. there may be one more 95-GHz 
component). 

The remaining spot B, is exceptional in the sense that all ob- 
served transitions have been detected towards it. Corresponding 
spectra are shown in Fig. \^ As in Fig. |5| the error bars show the 
3(T velocity uncertainty, with the exception of the high J 25-GHz 
transitions where the uncertainty is less than the spectral resolution, 
and therefore is not shown. The 95-GHz profile is wider than the 
84-GHz one and has an additional emission component. In contrast, 
all 25-GHz profiles are narrow. The 9.9- and 104-GHz profiles re- 
semble each other and have an almost single spectral channel spike 
on top of a broad symmetric line. 

To provide quantitative characteristics of the observed spectral 
profiles each of them has been fitted with a number of Gaussian 
components. The results of these fits are summarized in Table |2| 
As before, the uncertainties are given in parentheses and expressed 
in units of the least significant figure. The first column represents 
the name of the spot. The next six columns are self-explanatory 
and contain the parameters of each Gaussian component. The spec- 
tral and spatial domain were fitted separately. First, the peak radial 
velocity with respect to the local standard of rest (LSR), the line 
widths (FWHM) and the peak flux densities listed in the columns 2, 
5 and 6 were determined directly from the spectra shown in Fig.|2| 
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Figure 2. Spectra of spots A, C - F at 84 (upper plots) and 95 GHz (lower plots). Spots D and E were marginal detectons at 84 GHz. The error bar near the 
top of each plot shows the 3cr velocity uncertainty corresponding to the uncertainty in rest frequency (see TablelTl. 



and|3| Then, a 2D Gaussian source was fitted to the image of a 
single plane of the spectral cube, which corresponded to the peak 
velocity of each Gaussian component. This second fit gave the po- 
sition and the deconvolved source size estimate shown in columns 
3, 4 and 7. The latter is obtained under the assumption that the syn- 
thesized beam is known exactly and the source is circular. In the 
size column of Table |2| two components are described as 'unre- 
solved' (the fitted size in these cases was slightly smaller than the 



beam, owing to the effect of noise). Marginal 84-GHz detections to- 
wards the spots D and E are hidden in the edge noise in the image. 
Therefore, no second stage fit has been done for them. However, 
the position of these 84-GHz features, if they are real, is within 4" 
(the box size used in forming the spectra from the image cube) of 
the position of the appropriate 95-GHz spot. A combination of a 
5-function and a Gaussian was used to fit the 9.9- and 104-GHz 
spectra, which contain a narrow spike. An upper limit equal to the 
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Figure 3. Spectra of spot B, where masers in all observed transitions were detected. The error bar at the top of some plots shows the 3cr velocity uncertainty 
con'esponding to the uncertainty in rest frequency (see TablelTl. For J = 5 — 9 transitions at 25 GHz this uncertainty is below the spectral resolution and, 
therefore, is not shown. 



spectral resolution is given in the line width column of Table|2|for 
such narrow spikes. The flux density of these spikes exceeds the 
rms noise 8 and 22 times for the 9.9 and 104 GHz profiles respec- 
tively (see Table|2j- This fact leaves no doubts that the spikes are 
real. 

All uncertainties listed in columns 3-7 are the formal uncer- 
tainties of the Gaussian fit (lu). The formal uncertainty of the peak 
velocity is always half of that for the line width. However, due 
to blending of different components seen in most cases, there are 
likely to be systematic errors. Due to these errors, the individual 
components are likely to appear closer to each other according to 



the fit (given in Table|2j than they are in reality, both in the spatial 
domain and velocity. Their line widths and sizes are likely to be 
overestimated. In some cases, a superposition of two components 
is completely dominated by the brightest feature located at approx- 
imately the same velocity as the weaker one. We were unable to 
produce a reliable image plane fit for a weaker line in such situa- 
tions. 

The rest of the columns describe each line profile as a whole, 
providing more robust values than the results of the Gaussian fit. 
The columns are the peak LSR velocity and corresponding flux 
density, the area under the line profile (the integrated flux density), 
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Table 2. Fit results and profile parameters. The uncertainties are given in parentheses and expressed in units of the least significant figure. Notes: (a) the 
uncertainty is half of that for the line FWHM, (6) the uncertainty is the spectral resolution listed in TablelTI (c) the uncertainty is the rms listed in Tabled (d) 
the position is the same as for the component above as they can not be distinguished in the image, (e) the given upper limit is the spectral resolution, rather 
than a FWHM estimate. 
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52 28.68 (8) 0.28 (4) 


2.1 (3) 


0.9 


(2) 










B 






9-1 - 8-2 E (9.9 GHz) 








-31.56 


9.5 


2.80 (4) 


5.3 X lO'^ 




-31.572 


16.460 (2) 


52 25.73 (3) 0.32 (1) 


8.1 (3) 


0.10 


(9) 












-31.554'' 




<0.029'= 


2.4 (3) 




















22 - 2i E (25 GHz) 








-31.71 


0.3 


0.032 (3) 


2.8 X 10'* 




-31.69 


16.44 (2) 


52 26.4 (2) 0.15 (8) 


0.2 (1) 


3 


(2) 
















32 - 3i E (25 GHz) 








-31.61 


5.9 


1.303 (8) 


2.6 X 10^ 




-31.620 


16.459 (2) 


52 25.90 (4) 0.186(7) 


5.2 (2) 


0.08 


(3) 












-31.63'' 




0.46 (16) 


0.6 (2) 




















42 - 4i E (25 GHz) 








-31.71 


8.7 


2.17 (6) 


3.9 X 10® 




-31.715 


16.459 (2) 


52 25.90 (4) 0.203 (7) 


8.3 (3) 


0.24 


(6) 












-31.8'' 




0.5 (3) 


0.6 (3) 




















52 - 5i E (25 GHz) 








-31.69 


6.7 


1.80 (2) 


3.0 X 10® 




-31.678 


16.461 (8) 


52 25.70(11) 0.211 (5) 


6.4 (2) 


0.093 (8) 












-31.82'* 




0.3 (1) 


0.5 (2) 
















-31.24 


16.581 (6) 


52 25.5 (1) 0.5 (2) 


0.4 (2) 


0.5 


(4) 
















62 - 6i E (25 GHz), May 








-31.69 


7.5 


2.047 (4) 


3.3 X 10® 




-31.701 


16.464 (3) 


52 25.76 (5) 0.211 (5) 


7.3 (2) 


0.17 


(V) 












-31.89'' 




0.18 (6) 


0.5 (2) 
















-31.23 


16.594 (1) 


52 25.32 (2) 0.5 (1) 


0.6 (2) 


0.4 


(1) 
















62 - 61 E (25 GHz), June 








-31.67 


7.2 


1.88 (6) 


6.4 X 10® 




-31.649 


16.454 (1) 


52 25.70 (2) 0.203 (9) 


6.8 (3) 


0.079 (7) 












-31.80'' 




0.22 (9) 


0.8 (3) 
















-31.2 


16.581 (1) 


52 25.25 (1) 0.5 (3) 


0.5 (3) 


0.17 


(5) 
















72 - 7i E (25 GHz) 








-31.64 


6.8 


1.77 (4) 


5.9 X 10® 




-31.660 


16.462 (1) 


52 25.73 (1) 0.24 (1) 


6.6 (3) 


unresolved 












-31.3 


16.592 (2) 


52 25.26 (3) 0.3 (2) 


0.4 (3) 


0.28 


(5) 
















82 - 81 E (25 GHz) 








-31.67 


5.6 


1.48 (2) 


4.9 X 10® 




-31.644 


16.461 (3) 


52 25.78 (5) 0.24 (1) 


5.5 (3) 


0.060 


(8) 
















92 - 9i E (25 GHz) 








-31.63 


3.9 


0.956 (4) 


3.3 X 10® 




-31.62 


16.458 (1) 


52 25.70 (1) 0.24 (2) 


3.6(3) 


0.06 


(2) 
















5-1 -4oE(84GHz) 








-31.58 


12.3 


5.35 (2) 


3.7 X lO'' 




-31.85 


16.470 (8) 


52 26.29 (6) 0.38 (3) 


2.5 (2) 


0.8 


(4) 












-31.561 


16.454 (2) 


52 25.72 (2) 0.215 (5) 


9.3 (2) 


0.46 


(8) 












-31.31 


16.470 (5) 


52 25.73 (4) 0.59 (4) 


3.4 (2) 


0.68 


(8) 
















80 - 7i A+ (95 GHz) 








-31.52 


27.90 


19.0 (1) 


6.8 X 10'' 




-32.35 


16.469 (2) 


52 26.52 (2) 0.59 (6) 


2.8 (2) 


0.4 


(1) 












-31.808 


16.467 (5) 


52 26.15 (4) 0.415 (9) 


7.2 (2) 


0.7 


(1) 












-31.516 


16.463 (1) 


52 25.83 (1) 0.329 (3) 


20.3 (2) 


0.30 


(2) 












-31.235 


16.473 (3) 


52 25.86 (3) 0.594(9) 


10.4 (2) 


0.60 


(5) 












-30.56 


16.48 (1) 


52 26.1 (1) 0.46 (8) 


1.0(2) 


1.3 


(5) 
















11-1 - 10-2 E (104 GHz) 








-31.60 


9.59 


2.348 (3) 


2.0 X 10" 




-31.594 


16.462 (2) 


52 25.64 (2) <0.022'= 


2.2 (1) 


0.2 


(1) 












-31.580'' 




0.293 (5) 


7.5 (1) 















and the brightness temperature lower limit. The latter is calculated 
using the peak (in the spectral domain) of the total brightness con- 
tained in the same region of image, which was used to compute 
the spectra. This approach is more appropriate for unresolved or 
barely resolved sources and does not depend on (sometimes very 
crude) size estimates. Such high values of the brightness tempera- 



ture (lower limits from 6 x 10^ to 5 x 10*^ K, see Table|2j are gen- 
erally considered as an indication that the detections have a maser 
rather than thermal origin. The uncertainty of the peak velocity is 
equal to the spectral resolution, which is different in each particu- 
lar case (see Table0. A systematic offset of the velocity scale for a 
transition can be caused by rest frequency uncertainties, which are 
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Table 2 - continued 









Cjsussian components 








Peak 


Peak 








Spot 


LSR 


"2000 


(52000 Line 


Flux 


Size 


LSR 


flux 


//(-) 


dv 


n 




Velocity 


le'^ss™ 


—42 rWHM 


density 






velocity 


density 




Limit 




(kms-i) 




(arcmin arcsec) (kms~^) 


(Jy) 


(arcsec) 


/I — 1\ 
(km s ) 


(Jy) 


(Jy km s 




(K) 


C 


-31.14 
-30.77 


16.71 (1) 
16.672 (2) 


5-1 -4oE(84GHz) 
52 20.7 (1) 0.30 (5) 
52 20.62 (2) 0.28 (4) 


1.2 (2) 

1.3 (2) 


1.0 
0.4 


(1) 
(2) 


-30.77 


2.1 


0.7 


(1) 


6.0 X 10^ 








8o - 7i A+ (95 GHz) 








-30.76 


7.92 


3.549 (3) 


1.9 X 10"' 




-31.46 


16.75 (1) 


52 20.7 (1) 0.32 (5) 


0.7 (1) 


2.0 


(5) 














-31.167 


16.707 (2) 


52 20.85 (1) 0.270(8) 


3.5 (1) 


0.36 (5) 














-30.747 


16.686 (2) 


52 20.79 (1) 0.276 (4) 


7.8 (1) 


0.36 (3) 












D 


(marginal detection) 


5-1 -4oE(84GHz) 








-30.77 


3.0 


0.8 


(2) 






-30.8 




0.6 (2) 
8o - 7i A+ (95 GHz) 


1.5 (5) 






-30.51 


9.5 


7.54 


(4) 


1.7 X lO"' 




-30.92 


16.82 (1) 


51 59.63 (9) 0.29 (3) 


5.3 (5) 


0.7 


(6) 














-30.45 


16.850 (7) 


51 59.98 (5) 0.62 (4) 


9.0 (5) 


0.7 


(3) 












E 


(marginal detection) 


5-1 - 4o E(84 GHz) 








-30.33 


1.5 


0.4 


(1) 






-32.18 




0.09 (5) 


0.5 (3) 


















-30.336 




0.05 (1) 


1.2 (3) 


















-30.2 




0.9 (7) 
8o - 7i A+ (95 GHz) 


0.3 (3) 






-30.41 


7.1 


3.29 


(9) 


1.6 X lO"* 




-33.28 


17.53 (1) 


52 05.14(9) 0.16 (3) 


1.1 (2) 


1.0 


(5) 














-31.40 


17.72 (1) 


52 03.80(8) 0.12 (4) 


0.5 (2) 


2 


(2) 














-31.16 


17.48 (2) 


52 03.9 (2) 0.19 (5) 


0.8 (2) 


1.4 


(3) 














-30.403 


17.549 (4) 


52 04.09 (3) 0.210(8) 


5.1 (2) 


0.62 (3) 














-30.185 


17.551 (2) 


52 04.30(1) 0.32 (1) 


5.2 (2) 


0.38 (7) 












F 






5_i -4oE(84GHz) 








-27.28 


20.4 


9.3 


(2) 


5.8 X 10'' 




-27.87 


16.2 (1) 


52 00 (1) 0.18 (5) 


1.7 (4) 


unresolved 














-27.436 


16.163 (7) 


52 01.02 (6) 0.37 (1) 


11.1 (4) 


0.8 


(3) 














-27.311 


16.187 (8) 


52 00.99 (6) 0.106(5) 


9.8 (4) 


0.8 


(2) 














-27.108 


16.146 (5) 


52 00.73 (4) 0.254 (9) 
So - 7i A+ (95 GHz) 


12.8 (4) 


0.6 


(3) 


-27.14 


36.6 


20.56 


(7) 


8.7 X lO'' 




-27.94 


16.190(7) 


52 01.14(6) 0.10 (2) 


2.0(3) 


0.9 


(4) 














-27.75 


16.149(4) 


52 01.11 (3) 0.16 (2) 


3.5 (3) 


0.5 


(2) 














-27.321 


16.171 (3) 


52 01.23 (2) 0.482 (5) 


30.3 (3) 


0.53 (5) 














-27.094 


16.156(1) 


52 01.13 (1) 0.238 (5) 


16.9 (3) 


0.36 (4) 













also listed in Table Q This also applies to the peak velocity of the 
Gaussian components. The uncertainty of the integrated flux den- 
sity was obtained from either the noise level in the spectrum (see 
below) or the offset from zero of the mean flux density in the line- 
free part of the spectrum; we have quoted the larger of these two 
uncertainties in each case. 

To quantify the noise and provide useful upper limits for tran- 
sitions which have not been detected we have calculated the rms 
noise in the spectra constructed for each spot. The spectra for non- 
detections were extracted from the cube in essentially the same way 
as the spectra with a detected line emission shown in Fig. |2| and 
|3| For the spectra with a detected emission, rms calculations were 
confined to the emission free part of the spectrum. The la noise 
levels are listed in Table|3| They vary from spot to spot for the high 
frequency transitions due to the small primary beam size (compa- 
rable with the spread of the spots in the source) and the lack of 
an adequate mosaicing in the North-South direction. The noise lev- 
els list in Table|3|are the uncertainties of the peak flux densities in 
Table|2|(8th column). 

3.3 Notes on temporal variability 

Methanol masers are known to be vari- 
able (e.g. ICaswell. Vaile & Ellingsei] ll995iJ: 



iGoedhart. Gavlard & van der Walt! |2004) . although the sub- 
ject is not very well studied for class I masers. One of the goals 
of this work is to provide data for further development of maser 
models. The models are usually checked by comparing the flux 
density ratios of different transitions. However, these ratios can 
be affected not only by the studied pumping mechanism, but also 
by any temporal variability, if observations are not simultaneous. 
Temporal variability is typically not included in modelling and 
requires a long time series of data to be well understood. To obtain 
some estimate of the scale of temporal variability in the source we 
repeated the observations of some transitions in different observing 
runs, in two cases (see Table without a proper imaging (a 
single short scan). Provided the position of the source is known 
from another observing run, such non-imaging data still allows 
us to estimate the peak flux density and, therefore, to assess the 
variability. We have performed calibration similar to that used for 
the main dataset and obtained the following estimates for the peak 
flux density by fitting a point source at the known position to the 
visibility data: 6.7±0.2 Jy for the 25-GHz J=6 data taken in the 
August session, and 9.0±0.2 Jy for the 9.9-GHz data taken in May. 
The 25-GHz J=6 transition was observed in all three sessions, and 
the above value of the flux density is in agreement with that in 
Table |2| within 3a uncertainty. The 9.9-GHz non-imaging data are 
in agreement with the corresponding imaging experiment as well. 
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Table 3. The noise level (l-cr-rms) in the spectra. 





Transition 






Maser spot 








and 


A 


B 


C D 


E 


F 




frequency 






l-tr rms (Jy) 






9- 


i-8_2E (9.9 GHz) 






0.3 for all spots 






J2 


- Jl E, May (25 GHz) 






0. 1 for all spots 






J2 


- Ji E, June (25 GHz) 






0.2 for all spots 






5- 


1 - 4o E (84 GHz) 


0.2 


0.2 


0.2 0.5 


0.3 


0.5 


8o 


- 7i A+ (95 GHz) 


0.09 


0.09 


0.07 0.3 


0.2 


0.4 


11 


-1 - 10-2 E (104 GHz) 


0.1 


0.07 


0.08 0.4 


0.3 


0.5 



The latter is also important for our conclusion that polarization 
effects do not significantly affect the results of our observations 
(section |2j. This is because, in contrast to the June observations, 
the 9.9-GHz non-imaging observations carried out in May used a 
dual polarization mode and provided effectively a Stokes-I mea- 
surement. Taking into account the systematic errors in establishing 
the absolute flux density scale, which are estimated to be around 
10% and 3% at 25 and 9.9 GHz, respectively, we can conclude 
any temporal variability within the period of our observations is 
less than the calibration errors. Therefore, we can consider that we 
have inputs for the modelling which are free from any significant 
variability bias. 



4 DISCUSSION 

4.1 Association with outflow 

All the detected maser spots are located within a de nse molecular 
core (deconvolved FWHM 27", lOarav et alj|2003h and have ve- 
locities close to that of the ambient molecular material in the core 
(—30.0 km s~^). A remarkable correlation between the southern 
maser spots (labelled A to C) and the shocked gas mentioned in 
section im (see also Fig.0 points to an association of the class I 
masers with the interface region of this outflow. This result is some- 
what simi lar to that obtained for the n orthern star-forming region 
DR 21 by Plambeck & Menten ( 1990) more than a decade ago. In 
addition, Kurtz et al. ( 2004) revealed the same association for a few 
other northern sources. However a more favourable geometry and 
probably lesser complexity in the IRAS 16457-4247 star-forming 
region allows us to study this association in detail. 

In contrast to the southern spots, no associated H2 emission 
has been detected in the vicinity of the northern spots E and F 
(see Fig.0. However, spot D resides close to the northern con- 
tinuum lobe, which is associated with faint H2 emission. This 
morphology leaves no doubts that spot D is very similar to the 
southern maser spots, because the continuum source represents a 
working surface of the jet (.Brooks et al-.2003.,) . Moreover, it is 
evident from the large-scale image of the H2 emission (Fig. 
that an outflow lobe extending northwards is present. Therefore it 
is likely that all maser spots in this source are associated some- 
how with the molecular outflow. Two factors or their combina- 
tion can explain the non-detection of the H2 emission in the vicin- 
ity of the maser spots E and F. First, the 2.12 fim H2 emission 
can be affected by extinction, which may be higher in the north- 
em part of the region. Such var iations of ext inction hav e been 
observed in other sources (e.g. lEisloffel. Sm ith & Da vl3l2OO0l: 
iNanda Kumar. Bachiller & Davi j 120021) . In addition. iGara^^U 
(2006) have recently discovered a bipolar molecular outflow to- 
wards IRAS 16457-4247. This outflow is inclined to the line of 



sight at an angle of 40° with the receding flow located in the north- 
em part of the region. A higher extinction to the north from the 
central object has a natural explanation in this geometry. 

Alternatively, or in addition to this, the conditions in the inter- 
face region may be less energetic in the northern part of the source, 
if for example most of the dense material blocking the path of the 
outflow has already been swept away. The cooling timescale of the 
gas radiating shocked H2 e mission is very short and may be of the 
order of several years (e.g., 'Davis & SmitlJ ll995l : ISmith & Bran^ 
199C). This estimate is in agreement with measurements of vari- 
ability performed by Micono et al. ( 1998) for the HH 46/47 out- 
flow system. Because the IRAS 16457-4247 outflow extends well 
beyond the region containing maser spots (see Fig.0, the shocked 
gas in the vicinity of northem maser spots generated when the jet 
carved out a cavity in the dense molecular core, has already had 
enough time to radiate the energy and cool down if the interaction 
became less energetic for some reason. It is reasonable to expect 
that the shocks associated with intemal working surfaces (where a 
fast moving gas impacts a slow moving gas and drives shocks into 
the medium) are less energetic than those associated with the termi- 
nal working surfaces (where the flow impacts the cloud directly). 
To the south from the central source, the jet may have impacted the 
walls of the cavity generating a young curved shock. However, as it 
was mentioned above, there are a number of other factors that can 
account for the shape of this H2 feature. It must be noted that this 
scenario is speculative as there is no direct observational evidence 
at this stage that suggests a weaker interaction of the northern flow 
with the cloud. 

The main role of shocks for methanol maser formation is 
thought to be that they increase the methanol abundance in the gas 
phase. The production of methanol by purely gas-phase chemical 
reactions is inefficient, but the passa ge of weak shocks is a ble to re- 
lease it from the grain mantles (e.g. iHartauist et alJl995h . Shocks 
generated due to the interaction between outflows and the ambi- 
ent material are known to enhance the methanol abundance (e.g., 
Gibb & Davis 1998; Garay et al. 2002). In addition, outflows com- 
press the medium, which makes collisions more frequent and the 
pumping mechanism more efficient. The shocks have to be rel- 
atively gentle as methanol can survive sputtering or desorption 
of grain mantles only at shock velocities not greatly exceeding 
10 km s^^ (Garav et al. 2002). Therefore, class I methanol masers 
are likely to reside in the wakes associated with the bow-shocks or 
internal working surfaces rather than near the apex of a terminal 
bow-shock. The fact that all maser spots are only found close to the 
jet (Fig.0 supports this idea. 

It is worth noting that spot B, which is active in all observed 
transitions, is located near the edge of the brightest H2 knot. It is 
the only spot in this source where maser emission at 9.9, 25 and 
104 GHz has been d etected. According to models calculated by 
ISobolev et alj i2005h these masers require more energetic condi- 
tions (higher densities and temperatures) than the common masers 
at 84 or 95 GHz. Our results confirm this model, although with 
the caveat that the observed H2 emission may be affected by vari- 
able extinction. Spectroscopy of a number of H2 transitions, as well 
as observations of other shock tracers less prone to extinction (e.g. 
SiO and thermal methanol transitions) are required before firm con- 
clusions can be drawn. The 9.9- and 104-GHz profiles contain a 
narrow spike (Fig. |3j which has a width comparable to the spec- 
tral resolution (<0.03 km s~^, see Table0, and is the narrowest 
spectral feature ever claimed according to our knowledge. Such a 
narrow feature implies that the maser is unsaturated and requires 
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Figure 4. The correlation between flux densities at 84 and 95 GHz. Each 
data point is labelled with the spot name and velocities in km s~^ (two 
numbers in brackets) corresponding to the 84- and 95-GHz components 
taken from Tablel2l The error bars correspond to the 3a" formal uncertainty 
of the Gaussian fit. The dashed line with the slope 2.1±0.3 represents the 
least square fit discussed in the text. 



Fig. |4] shows the peak flux densities {F) at 84 and 95 GHz 
for matching pairs of the Gaussian components listed in Table |2| 
Each point is labelled with the spot name and the peak velocities of 
the corresponding 84- and 95-GHz Gaussian components. The er- 
ror bars show 3(t formal uncertainties of the fit and do not take into 
account possible systematic errors due to the absolute flux density 
scale and the structure of each spectral profile. It was concluded 
from the general appearance of the profiles that the spot F is likely 
to have an additional component at 95 GHz. Therefore, the system- 
atic errors are particularly acute for this spot. Despite this caveat, 
Fig.|4|reveals the correlation between flux densities quite well. A 
linear least square fit gives the following dependence 

F(95 GHz) = (2.1 ± 0.3) x F(84 GHz) (2.6 ± 2.5), (1) 

which is shown in Fig. |4] by the dashed line. The additive term 
is not statistically significant. The absolute flux density scale cal- 
ibration affects all spots at the same frequency in the same way. 
The refore, the slope has an additional uncertainty of 30% (about 
0.6'). lKalenskii et ai]i200ll) obtained a similar result to Q by com- 
paring single dish surveys of a large number of sources at these 
two frequencies. However, to study such a dependence reliably, an 
interferometer is required, due to the spread of maser spots. It is 
worth mentioning that interferometr ic observations of the D R21 
star-forming complex carried out by IB atria & Mentenl il988 h and 
f'lambeck & Menten ( 1990), found the ratio of the 95- to 84-GHz 
flux density to be approximately 5.5, somewhat higher than ob- 
served in IRAS 16547-4247. 



turbulent motions to be suppressed (e.g. by the magnetic field) in 
the gas involved in the maser action. 

The clear association of class I methanol masers with outflows 
in this and other sources suggests that searches for class I methanol 
masers should be extended to sites of intermedi ate- and low-mass 
star f ormation, where outflows are common iReipurth & Ballvl 
I2OOII) . Indeed, at such si tes a few class I masers have recently been 
found tealenskii et'arj|20ol) . but a lack of accurate positions for 
these masers means their association with outflows cannot be veri- 
fied at the present. 

4.2 Masers at 84 and 95 GHz 

As mentioned above, the maser profiles at 84 and 95 GHz look 
similar in general (Fig.|2|and|3}- This suggests a com mon range of 
physic al conditions required for these masers to form. lVal'tts et alJ 
yOOy) found a similar correspondence between the 95-GHz masers 
and the 7o — 61 A+ masers at 44 GHz (this frequency is currently 
outside the band of the ATCA receivers and has not been observed 
in this project) using a large sample of southern maser sources, 
including IRAS 16547-4247. The flux densities of these two tran- 
sitions were found to correlate, despite a significant scatter of the 
points, which is most likely caused by the lack of accurate positions 
for the majority of known class I methanol masers. However, the 
44- and 95-GHz masers belong to the same (J + l)o — Ji A^ series 
of transitions (J=6 and 7 respectively) and are expected to correlate. 
In contrast, the 84 and 95-GHz maser transitions belong to a differ- 
ent series of transitions, and even to different species of methanol 
(E- and A-methanol respectively, see Table 0, which act as two 
different molecules. Therefore, a correlation between the flux den- 
sities corresponding to these two transitions is not self-evident and 
can be a constraining factor for models of maser pumping. 



4.3 Rest frequencies 

High spectral resolution observations can be used to improve the 
accuracy to which the transition rest frequencies are known. How- 
ever, in many cases, the laboratory measurement is preferred re- 
gardless of its accuracy because the assumption of equal velocities 
for all transitions may be misleading (e.g., the discussion about the 
95- and 25-GHz masers in the next section). According to Fig. |3| 
the 9.9- and 104-GHz transitions have a single channel spike on 
top of the broad line. These two transitions belong to the same 
(J-l-l)_i — J_2 E transition series (J=8 and 10) and are expected 
to show a similar behavior Therefore, assuming that the spike is at 
the same velocity, the rest frequency of the 9-1—8-2 E transition is 
9936.201 ±0.001 MHz. The uncertainty is determined by the spec- 
tral resolution and the rest frequency uncertainties of the 9.9- and 
104-GHz transitions. The rest frequency was adjusted with respect 
to TableQby 1 kHz or 0.03 km s~^. It is not practical to adjust the 
rest frequency for other transitions. 

4.4 Maser models 

Theoretical studies of class I methanol masers were impeded for 
a long time by the absence of an adequate model of collisions for 
the methanol mo lecule. Recently calculated rate coefficients (e.g., 
IPottage. Flower & Davis 2004) ha ve enc ouraged pumping stud- 
ies of these masers. ..Sobolev et aP i2005l) performed a qualitative 
analysis and revealed that the 25-GHz masers are bright in mod- 
els with specific column density (equivalent to the column den- 
sity divided by the line-width) greater than 10^^ cm^^ s and re- 
quire relatively high kinetic temperatures (75 — 100 K) and densi- 
ties (10^ — 10^ cm~^). Similar conditions are required for the 9.9- 
and 104-GHz masers, although they prefer a greater beaming (elon- 
gated geometry) and either lower densities or higher temperatures 
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Table 4. Model results for spot B. ij is the ratio of the integrated flux density to that of the J=4 transition at 25 GHz, t is the optical depth and Tj, is the 
brightness temperature. The uncertainties of the observed log(ij) are given in the parentheses and expressed in units of the least significant figure. According 
to Tablel2lthe brightness temperature of the J=4 25-GHz transition exceeds 3.9 X 10'' K. 



Transition u Observed Model 1 Model 2 

(GHz) log(i?) log(T6,K) r log(ij) log(ri,,K) r log{i?) 
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(more than 100 K). Spot B is an excellent test site to develop a 
class 1 pumping model because of the large number of transitions 
detected. However, a detailed investigation of the vast parameter 
space of the maser models is beyond the scope of this paper and will 
be reported elsewhere. At present we are unable to reproduce the 
relative brightnesses of all transitions detecte d in spot B in a single 
model similar to that o f'Sobol eveta'i] j2005'). We suspect that gra- 
dients in the physical parameters may be influential and leading to 
the co-existence of the different pumping regimes on spatial scales 
unresolved by ATCA. It is expected that strong gradients may ac- 
company shocks, although their effects on the maser pumping are 
yet to be studied. According to Table|2|the peak components at 25 
and 95 GHz are displaced by not more than O'.'l. Most likely this 
displacement is caused by a systematic error due to the structure 
of the 95-GHz maser because no such displacement is observed for 
the 84-GHz maser. Taking this uncertainty as an upper limit and 
assuming a distance to the source of 2.9 kpc, the parameters must 
change on a scale around 300 AU or less. 

It is worth noting, however, that firm conclusions about any 
gradients can only be made after an extensive parameter search. 
In this paper we have just illustrated the situation using two repre- 
sentative models sel ected from the range of models calculated by 
ISobolevetal] j2005). Outlines of the modelling procedure applica- 
ble to class I masers ca n be found in^ronkov et al. t.2005b.) and 
the code is described bv lSutton et alJ^OCMl) and references therein. 
Both models use a kinetic temperature (Tgas) of 75 K, an intrinsic 
linewidth (to simulate blending of spectral lines) of 0.5 km s^^, 
and a methanol abundance relative to hydrogen of 10~^. The pa- 
rameters that differ between these two models are the hydrogen 
number density (n/j), 10^ and 10^ cm~"^, the specific column den- 
sity (equivalent to the column density divided by the line-width), 
10^^ and lO^''-^ cm"^ s, and the beaming factor (e^^), 20 and 
30, for the first and the second models respectively. The results 
of the simulations are shown in Table |4] Most of the columns are 
self explanatory, the columns labelled log R give a decimal loga- 
rithm of the ratio of the integrated flux density of the transition to 
that of the J=4 25-GHz transition (as observed and for each of the 
models). The uncertainties of the observed log(_R) are given in the 
parentheses and expressed in units of the least significant figure. 
These are cumulative uncertainties, which, in general, account for 
both the absolute flux density scale uncertainty and that of the inte- 



grated flux density given in Table|2| The J=2 and 3 transitions were 
observed simultaneously with the reference (J=4) transition in the 
same band (see section |2j. Therefore, the ratio R is not affected 
by the absolute flux density scale calibration for these two transi- 
tions. Other columns contain the brightness temperatures (log Tt) 
and optical depths (r) predicted by the two models. 

Figure|5|illustrates the dependence of log R on the excitation 
energy for the 25-GHz series (data from Table |4}- Open squares 
and triangles correspond to the first and the second models, respec- 
tively. Observed values are shown by the filled circles with error 
bars corresponding to the 3a uncertainty. The observed dependence 
represents essentially a rotational diagram (because the product of 
the spontaneous decay rate by the statistical weight does not vary 
much throughout the series), which is a well k nown method to ana - 
lyze the level populations of the molecule (e.g. lMenten et alJl98d) . 
The plot of observed values shows a significant curvature, although 
a linear fit to the highest J transitions is possible. Such a fit, which 
takes into account all transitions with J> 3, has a rotational temper- 
ature of 80±10 K. It should be remembered that this temperature 
is unlikely to be related to any physical or even excitation tempera- 
ture. The latter is estimated correctly if the fitted 25-GHz transitions 
are all optically thin or have the same optical depth. 

It is evident from Table |4] that despite the fact that the first 
model describes the ratios for the 9.9 and the 25-GHz transitions 
quite well, it fails to produce the masers at 84, 95 and 104 GHz. 
Similarly, the second model is adequate for 9.9-, 84-, 95- and 104- 
GHz masers, but fails to describe the behavior of the high J 25- 
GHz transitions (Fig.|5}- The brightness temperature of the refer- 
ence 42 - 4i E transition is 1.9x 10** and 3.6 x 10^ K, for the first 
and the second models respectively (see Table|4}. Both these val- 
ues are in agreement with the observed lower limit (3.9 x 10^ K, 
Table|2j. 



4.5 Evolutionary stages and maser activity 

H2O, OH and both classes of methanol masers often, but not al- 
ways, are present in the same region of high-mass star forma- 
tion. The question whether these masers trace an evolutionary se- 
quence (with a pa rtial overlap in time) has been in vestigated in var- 
ious studies_^e;2jForst^^^Cas & MoscadeUil 
|200(1 ISzvmczak. Pillai & MenterJ l2005t lEllingser< 120051 12003^ 
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Figure 5. The dependence of the observed and predicted relative integrated 
flux densities (R) versus the exitation energy for the 25-GHz masers de- 
tected in spot B. Open squares and triangles correspond to the first and the 
second models (see TablePl. respectively. Observed values are represented 
by the filled circles. The J=4 transition is taken as a reference (see text). 
Therefore, the observed point and the model points coincide exactly for this 
transition. The error bars for observed values represent a cumulative 3(T 
uncertainty, which takes into account the uncertainty of the absolute flux 
density scale for J>4 transitions (J=2,3 and 4 transitions were observed si- 
multaneously in the same band; see section|5). There are two independent 
observations of the J=6 transition. 



From their analys i s of a large sample of star-forming regions 
Iporster & Caswell il989t) concluded that the compact, isolated 
H2O masers mark the earliest evolutionary stages, and OH masers 
appear afterwards. The class II methanol masers at 6.7 GHz seem to 
largely ov erlap w ith the H2O masers and precede OH (e.g., Szym- 
czak et al. l2005h . The majority o f class II masers are not associ- 
ated w ith detectable Hll regions jPhillips et alJ ll998HWalsh et alJ 
[199^, but show an association wit h millimetre and submillimetre 
Tpestaiozzi et a l.' 2002. Walsh et al. 2003) sources which leaves no 
doubt that they trace a very early stage of the (proto)stellar evolu- 
tion. However, the question as to where the class I methanol masers 
fit into this picture is still poorly understood. The main obstacle is 
the lack of untargeted surveys of class I masers, which have not 
been undertaken mainly because there is no widespread, bright, 
low frequency maser transition similar to the 6.7 GHz class II tran- 
sition, or mainline OH. Most class I masers, except for a few fa- 
mous sources like OMC-I, have been found towards known class II 
methanol or OH ma ser sites, which are often located within a sin- 
gle dish beam (e.g.. (sivsh et al . 1994; EUingsen 2005). Therefore, 
all these samples are biased. .EUingsen (.2006.) investigated this 
question using the infrared properties of a subsample of methanol 
masers associated with GLIMPSE (Galactic Legacy Infrared Mid- 
Plane Survey Extraordinaire) catalogue point sources. Despite the 
small number statistics, which makes any firm conclusions impos- 
sible, it was found that the sources with a class I maser seem to have 
redde r GLIMPSE colo urs than those without it. Based on these re- 
sults, lEUingserl i20o3) concluded that the class I methanol masers 
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Figure 6. The same as the inset in Fig.0 but with the GLIMPSE 5.8-/jm 
emission shown by greyscale. 



may signpost an earlier stage of high-mass star formation than the 
clas s II masers. 

iForster & Caswell il989h detected both mainline OH and 
H2O masers in IRAS 16547-4247 at a position close to the central 
continuum source in Fig.Q Taking into account the information 
given in the previous paragraph, one may expect the source to be 
relatively old and most likely to be active in class II maser tran- 
sitions. It is generally assumed that each high-mass star-forming 
region has a detectable 6.7 GHz maser sometime during the early 
stages of its evolution regardless of the maser variability and the 
source geometry (e.g., van der Walt 2005). However, no 6.7-GHz 
emission has been d etected in t his sou rce (upper limit of 0.3 Jy, 
ICaswell et aljl995bl;lwalsh et aljl998h . in contrast to the detection 
of the large number of class I masers reported in this paper. This 
fact encourages us to revisit the question of where class I methanol 
masers fit into an evolutionary sequen ce. 

The GLIMPSE survey ( Benja min et aljl2003h provides mid- 
infrared images of the Galactic Plane with unprecedented sensitiv- 
ity. Figure|5|shows the locations of the radio continuum, methanol, 
OH and water masers superimposed on the GLIMPSE 5.8-/^m 
emission (greyscale) on the same scale as Fig. Q The most strik- 
ing feature of this image is that there are no MIR sources asso- 
ciated with either the jet source, nor the class I masers. Diffuse 
emission, probably arising from the Polycyclic Aromatic Hydro- 
carbons (PAH) band at 6.2 /im appears to be associated with the 
outflow region and traces a similar area to the shock-excited H2. 
Although, unlike the H2, the PAH emission is stronger closer to the 
central continuum source and very faint in the vicinity of spot B. 
The GLIMPSE 3.6-/im image shows similar diffuse emission, also 
likely to be associated with PAH emission, in this case from the 
3.3-/im band. Fig.|6|also shows a strong, extended MIR source ly- 
ing between spots D and F. We have also examined the MSX im- 
ages of this region which show that although the I2-/im emission 
(C-band) is centred on the strong MIR source seen in Fig.|6| the 
21-/im emission (E-band) is offset and centred on the central con- 
tinuum source. That the source driving the jet/outflow/masers is 
only visible at wavelengths longer than 12 /j m suggests that this 
source is at a very early evolutionary phase. iGarav et'al] i20o3) 
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argued that the observed cm-wavelength continuum emission to- 
wards IRAS 16547-4247 is too weak for such a luminous object 
(6.2 X 10* L(7)). Su pported by the characteristics of the spectral 
line profiles TGarav et al. ( 2003) explained this flux discrepancy by 
the ongoing intense accretion, which forbids the development of 
an ultra-compact Hll region. Intense accretion is further evidence 
for an early evolutionary phase. Ellinesen (2006) found the vast 
majority of class II methanol maser sources to be associated with 
MIR emission visible in the GLIMPSE observations. That 

such emission is absent in this source argues that it is likely to be 
at an earlier evolutionary phase than that which produces class II 
methanol masers; however, this is in direct contrast with the pres- 
ence of OH masers which would generally be taken to support the 
opposite interpretation. 

In theory, a moderately dense gas (dense enough to switch 
the collisional pumping on, but not dense enough to produce 
thermalization), rich in methanol and separated from strong in- 
frared sources is able to produce most class I methanol m asers 
jCragg et aljll99i IVoronko^ll999t Rtonkov et alj|2605bl) . The 
fact that the class I masers do not require infrared radiation for 
pumping suggests that they might exist at earlier evolutionary 
stages than the class II masers, which need a nearby infrared source 
for pumping. However, a mechanism able to increase the methanol 
abundance above that seen in cold cores must already be active or 
have been active in the past to produce a detectable maser. As was 
discussed above, a release of methanol from the grain mantles by 
weak shocks is a good example of such a mechanism. One cannot 
exclude the possibility that class I masers survive even when the 
shock, which had increased the methanol abundance in the past, 
has already dissipated. Rare class I masers at 9.9 and 104 GHz 
require more energetic conditions jSobolev et aT]l2005h . They are 
most likely (or even exclusively) located near the young shocks. 
Spot B in Fig. Q is likely to be a good example of such a situ- 
ation. If the idea that the common class I masers can outlast the 
shock is correct, these masers become hard to destroy, once they ap- 
pear. Situated at a considerable distance from the central object they 
probably fade out slowly as the methanol abundance decre ases (e.g. 
due to freezeout or gas-phase chemical processing, Averv & Chiad 
Il996h . unless the chemical composition of the medium is altered 
once again by new energetic outflow events or by a slight change 
in the outflow direction. A new epoch of interaction between the 
outflow and the medium can either destroy the molecules engaged 
in the maser emission or produce new masers. Assumi ng typical 
outflo w velocities of 10^ — lO'' km s"^ (e.g.. Rei purth & Ballvl 
1200 ll) . it would take at least lO^-IO^ years for this disturbance to 
reach the maser sites for the given distance to the source (2.9 kpc) 
and the angular separations of masers from the central object (see 
Table |2j. A timescale for methanol destruction by ion-molecule 
reactions has approximately the same order of magnitude (e.g., 
ISandford & Allamandoliil993llGarav et all2002h . 

The class II masers are usually located very close to the YSO 
and, therefore, can be more easily destroyed during the course of 
its evolution. Given the detection of OH masers, which suggests an 
evolved stage, it is, therefore, possible that the 6.7-GHz maser has 
already been quenched in the studied source. We cannot presently 
resolve the contradiction between the infrared data and the presence 
of OH masers and confidently determine whether the source is too 
young or too old to have a 6.7-GHz maser. However, these spec- 
ulations show that the evolutionary stage where the class I masers 
are active, may in fact last longer and start earlier than when the 
class II masers are active. 



5 CONCLUSIONS 

(i) The class I methanol maser emission consists of a cluster of 6 
spots spread over an area of 30" in extent. Five spots were detected 
in only the 84- and 95-GHz transitions (for two spots the 84-GHz 
detection is marginal), while the sixth spot shows activity in all 12 
observed transitions. 

(ii) The three most southern maser spots show clear association 
with a jet-driven molecular outflow. Their velocities are close to 
that of the molecular core within which the jet is embedded. This 
fact supports the idea that the class I masers reside in the interface 
regions of outflows. 

(iii) Comparison with OH masers and infrared data reveals a po- 
tential discrepancy in the expected evolutionary state: the presence 
of the OH masers usually means that the source is evolved, but 
the infrared data suggest otherwise. The class II methanol masers 
at 6.7 GHz are very widespread and are usually treated as an in- 
evitable phenomenon at some stage of high-mass protostellar evo- 
lution (see section |431 . The lack of such a maser in this source 
(which is active in many class I maser transitions) raises an ad- 
ditional question concerning its evolutionary state: whether this 
source is too young or too old to have a 6.7 GHz maser? We ar- 
gue that both cases are possible because the class I masers do 
not require an infrared source for pumping and are harder to de- 
stroy than the class II masers, which are usually located close to 
the YSO. We, therefore, suggest that the evolutionary stage where 
the class I masers are present may last longer than that with the 
class II masers, although the exact status of this source is not clear 
at present. 

(iv) We report the first interferometric observations of the rare 
9.9- and 104-GHz masers. The spectra for these transitions contain 
a very narrow spike (< 0.03 km s~^) which has a brightness tem- 
perature greater than 5.3x10^ and 2.0 x lO"* K at 9.9 and 104 GHz, 
respectively. 

(v) High spectral resolution data leads us to suggest that the 
rest frequency for the 9_i — 8_2 E transition should be refined 
to 9936.201±0.001. 
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